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FURTHER NOTE ON MACLAURIN'S SPHEROID* 

BY 

GEORGE HOWARD DARWIN 

In volume 4 (1903) of these Transactions a paper of mine was published 
on the approximate determination of Maclaurin's ellipsoid. It was there shown 
that it is possible to determine the form of the ellipsoid and its angular velocity 
by means of spherical harmonic analysis to a higher order of approximation 
than is usually supposed to be attainable by that method. 

It appears that we ought to be able by means of the considerations adduced 
to determine the angular velocity corresponding to a given ellipticity as far as 
the cube of that ellipticity, and I was much puzzled to explain why I failed to 
obtain correctly the term involving the cube. 

On page 117 of the paper I gave reasons for believing that the determina- 
tion of this last term would prove to be illusory. Having had occasion recently 
to examine carefully another piece of work I saw that the reason given was in 
its turn fallacious, and that this last term should be attainable. Accordingly I 
looked carefully over the analysis to see whether there were any terms missing 
in my calculation, and having found the source of the discrepancy I propose in 
this note to supply the defect. No attempt will be made to make this new 
investigation complete, so that the present paper will only be intelligible in con- 
junction with the previous one. 

The equation to the spheroid which is to be made a figure of equilibrium 
under rotation a> was written 

T=-e£ 2 -/£,-£/• #', 

where t = (a 3 — ^)/3a 3 , S 2 = J — p 2 , and S t and S' t are the spherical surface 
harmonics of fi the cosine of colatitude from the axis of rotation. 

On page 118 I showed that the potential of the sphere from which the ellip- 
soid is derived is 

and I put r 2 = a 2 ( 1 — 2t — t 2 ) . Now one more term in the development of 
r 2 should have been included, and I ought to have written 
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r- 2 = a 2 (l-2T-T 2 -|T 3 -..). 

In order to evaluate the additional term in the part of the energy denoted 
— SH another integral is required, besides those given on page 117. 
It is, say, 

By the results given on page 117 for ( S 2 ) 2 we have 

_K) 2 , K) 2 K) 2 , 

4>o 4>z <t>i 

With the numerical values given on the same page, I find 

16 

a *~¥J7f 

Returning to page 118 it is easy to see that the additional term in SS is 

Hence the additional term in equation (5) which gives the value of %SS — SJR 
is(Jf 2 /a)(-*eVJ. 

On page 119 I avowedly dropped a term in CM, viz.: 



9M* 

iira 



Next on page 127 I avowedly dropped a term in the energy \DD, viz.: 

Lastly on page 129 I omitted the whole of the energy given by the formula 
(12). It is 

ipr 9 k) 2 5KH 

4 20 <f> 2 6 h J 6 ' 
Thus all the omitted portions together are 

a [3.5 °* 20 <£ 2 6 £ J 
With the numerical values given above the additional term in the lost energy 
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of the system (omitted because I thought terms in e 4 would lead to no further 
accuracy) is found to be 

a I 3\V e J' 
In calculating the moment of inertia C r on page 129, I omitted purposely 

Hence this with its sign changed and multiplied by Ja> 2 is an additional term in 
(14) giving jCtD 2 . Thus (14) should run 

It is convenient to insert numerical values throughout, instead of retaining 
symbols as I did before. 

I find then that the first formula of § 7 will run 

F .M*_ 4 _64_ 16 104 1 f-1 

.a^. --3 >5 2« -335276 -337726 +325276/ -32/ -i 2J+1 UJ P< 



o 



w 



n 



X7 e3 + O e ^]- 



+ 47rp[5 + 5 C + 3.7 e ~3 3 

In this the terms in e 4 in the first part and in e 5 in the second part are the 
additional terms now included. 

Since d_E/df. = is one of the conditions for the figure of equilibrium, it 
follows at once that_/"j = 0, for all values of i and s. 

By neglecting all terms above those of the second order, and putting 
dEjde = 0, we see that as a first approximation aP/iTrp = 4e/3.5 On putting 
dEjdf= we obtain 

104 2 2 tf_ _8 

3 a ;5 2 .7 e ~ & f+ 4tt P 5.7 e 

Then by means of the first approximation for &> 2 /47r/>, we find 

On substituting this value for f in the expression for the energy, we have 

¥1^ A 2 6 1_ 3 „ 35 ^„ 4 j^. T 3 2 ■ A 2 808 _ si 

■" "*" „ OC2 e OS C2 T e 03C2-r2 e +/( I E + C e + 3 J e + Q3 C >2 e - 



3.5 2 3 3 .5 2 .7 3 3 .5 2 .7 2 T 47rp [ 5 T 5 



*In line 7 from foot of page 129 there is a misprint ; the double integral should obviously be 
a single one, and dr should be deleted. 
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On equating d_E jde to zero we find 

/ _S 176 2 \_^f T 20_ 404 1 

V + 3.7 e + 3 2 .7 2e ) ~ i-rrp [ + 3.7 6 + 3 2 .7 2<5 J " 



K5 e{ 



4 

Whence 

-"--* e(l-*e + ^eA 
47rp~3.5 < V 7 e + 3.7 2 e )' 

Only the first two terms were evaluated in the former paper, and this is the 
required new result. 

The equation to the ellipsoid, as calculated in this way, is 

£=1+8(6$+^). 

The method then will give the coefficient of S i only as far as e 2 , and fails to 
give the coefficient of S e . Nevertheless it gives the square of the angular 
velocity as far as e 3 . 

It remains to verify the correctness of our result. 

It may be proved, as in the former paper, that the equation to an ellipsoid 
whose equatorial and polar semi-axes are a { and a x ( 1 — e l ) may be written in 
the form 

r s =a?(l-«i)[l+3(«, + ^?+A«?)^ + ¥«?(l+T 9 l«.)^-lK^]- 
Now a denotes the mean radius of our spheroid, so that 

a 3 = a\ ( 1 — e, ) . 

Since 3e, 3/", and say Sg are to be the coefficients S 2 , S 4 , S 6 , in the equation 
to the ellipsoid 

r 3 = a 3 [ 1 + 3e# 2 + 3/£ 4 + 3^] , 
we have 
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By 


inverting 


the 


expression 


for 


e, we find 
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If the eccentricity of the ellipsoid be denoted by sin 7 , we have 
■7=1— "n sm 2 j = 2e l (l—^e l J. 

; 7 sin 2 7 = 2e, M - g e, + g eM = 2e M - y e + -^j 2 e 2 J , 



COS 1 

Hence 



COS' 

/ 19 \ 

cos 7 sin 4 7 = 4e|( 1 — 2eJ = 4e 2 ( 1 — -y e I , 

cos 7 sin 6 7 = 8el = 8e 3 . 
Now the rigorous solution for the angular velocity of Maclaurin's ellipsoid is 
{In- l)!sin 2 "7 



= cos 7 ^ 



4tt P ~ ' V(2n+l)(2w + 3)[(n-l)!] 2 2 2 



2 3 5 

= g-g cos 7 sin 2 7 -f =-= cos 7 sin 4 7 + 02 yf cos 7 sin 6 7 



Substituting for 7 in terms of e , we find 

_ _4 / 4 ±_ \ 



4:TTp 

This agrees, as it shoidd, with the result found by means of spherical harmonic 
analysis. 

I do not think that the additional term found suffices to assign a limit to the 
stability of the ellipsoid, as I conjectured that it would do. 
Cambridge, England. 



